Local stabilizer codes in three dimensions without string logical operators 



Jeongwan Haatj^] 

Institute for Quantum Information, California Institute of Technology, Pasadena, CA 91125, USA 

(Dated: 28 February 2011) 

We suggest concrete models for self-correcting quantum memory by reporting examples of local 
stabilizer codes in 3D that have no string logical operators. Previously known local stabilizer codes 
in 3D all have string-like logical operators, which make the codes non-self-correcting. We introduce 
a notion of "logical string segments" to avoid difficulties in defining one dimensional objects in 
discrete lattices. We prove that every string-like logical operator of our code can be deformed to a 
disjoint union of short segments, each of which is in the stabilizer group. The code has surface-like 
logical operators whose partial implementation has unsatisfied stabilizers along its boundary. 
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I. INTRODUCTION 

Self-correcting quantum memory is an interesting sub- 
ject not only because of its application for quantum in- 
formation processing technology, but also because of its 
implication for quantum many-body physics; it shows 
a topological order at finite temperature. It is known 
that in 4D a self-correcting quantum memory is possible: 
Toric code [3 which is a CSS stabilizer code. There 
are classes of models in 2D that are not self-correcting 
[3H5] including those based on local stabilizer codes. It is 
thus a natural question whether a self-correcting quan- 
tum memory is possible in 3D, at least in the class of 
models based on stabilizer codes. 

A string-like logical operator plays an important role 
in the thermal instability; its existence is crucial in the 
no-go theorems [3H5] for self-correcting quantum memory 
in 2D based on local stabilizer codes, and more gener- 
ally, on local commuting projector codes. The string-like 
logical operator arises easily under the interaction with 
thermal bath, and hence, adversely affect encoded in- 
formation. Known models in 3D e.g., toric code OE], 
Chamon model [8^, topological color code [9" and Kim 
model [10 , do have string-like logical operators. Bacon 
subsystem code in 3D [TT] which does not have string-like 
bare logical operator, might be self-correcting, but it is 
not yet affirmative since its Hamiltonian is hard to solve. 

There is an issue of defining string-like logical oper- 
ator. Since a lattice is a discrete space, it is generally 
not possible to define the dimension of a subset of the 
lattice. An observation is that a string is a union of seg- 
ments, each of which has two end points. Thus we define 
logical string segments as a finite object that has two an- 
chors at the end with its middle part commuting with 
stabilizer generators. A string logical operator is then 
a logical operator that contains arbitrarily long logical 
string segments. 

The main result of this paper is that there exist local 
CSS stabilizer codes in 3D that are free of string logical 



* jwhaahOcaltech.ediil 



operators. We give the complete classification of codes 
{cubic codes) under our consideration in Sec. |ll] We ex- 
plain how we classify them in Sec. |III| In Sec. |IV[ we 
prove that the code distance is at least linear in system 
size. Sec. |V]is the central section where we define a log- 
ical string segment and prove that four of our codes are 
free of string logical operators. We report exact empiri- 
cal formulae of the number of logical qubits of our codes 
in Sec. |VI[ Finally, we discuss thermal stability of our 
codes and related issues in Sec. IVIII Sec. IVIIII contains 
our concluding remarks. 

Let us review the formalism of stabilizer codes. Let Vn 
be the group of Pauli operators acting on n qubits. An 
abelian subgroup S of Vn is called the stabilizer group 
if — / ^ S. The stabilizer group S defines a subspace of 
n-qubit Hilbert space by 

C = {IV') : sIV-) = |^)for ah s e S}, 

which is the code space. C is nonzero because — / ^ S. A 
CSS code is defined by a stabilizer group, each element 
of which can be written as a product of X- and Z- type 
stabilizer elements. The Pauli group has a nice property 
that any pair of elements is either commuting or anti- 
commuting, and that every element squares to identity. 
If we abelianize the Pauli group Vn by ignoring all phase 
factors |12| . we obtain 2n-dimensional vector space over 
the binary field equipped with a symplectic bilinear form 
A; A(a, 6) = 1 if a and b anti-commute, and A(a, 6) = 
if they commute. The product of two Pauli operators 
is expressed by the addition of the two corresponding 
vectors. The identity operator is the zero vector. 

In this respect, the stabilizer group is characterized 
as an isotropic subgroup. Note that the condition that 
—I^S should be checked separately. We abuse the no- 
tation and use the same symbol S to denote the vector 
space corresponding to the group. The orthogonal com- 
plement 5^ of S with respect to the symplectic form 
is the space of logical operators. The set of nontrivial 
logical operators modulo stabilizer group is the quotient 
space S-^/S. A stabilizer code is (geometrically-) ^oca^ if 
its stabilizer group is generated by (geometrically-)local 
Pauli operators. 
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A translation-invariant local stabilizer code can be de- 
fined on the infinite lattice. In this case, we define the 
stabilizer group as a group of all finite products of local 
generators. A logical operator is a Pauli operator possi- 
bly with infinite support that commutes with every gen- 
erator. Since each generator is local, the commutation 
relation between the generator and an arbitrary Pauli 
operator is well-defined. 



II. COMPLETE LIST OF CUBIC CODES 

We seek for a simple local stabilizer code that is 
translation-invariant, encodes at least one logical qubit, 
has large code distance, and does not have any string 
logical operator. (Formal definition of string logical op- 
erator will be given in Sec.|v]) To start with, consider a 
local stabilizer code on a D-dimensional simple cubic lat- 
tice with one qubit at each site. A general stabilizer 
generator may act on a bounded number of qubits in an 
arbitrary way. However, if we coarse-grain the lattice, or 
equivalently, put m > 1 qubits at each site of the lattice, 
we can say without loss of generality that a generator 
acts on the qubits on 2^ sites of a unit hypercube. 

We focus on stabilizer codes with only two types of 
generators for simplicity. Since each generator can be 
described by a 2m x 2^-component binary vector, there 
are a finite yet large number of conceivable generators. 
We will demand a certain structure of generators in order 
to reduce the number of candidates. We will further im- 
pose conditions such that the code does not contain any 
nontrivial logical operator on a "straight line" , which will 
be necessary for the codes to be without string logical op- 
erator. 

The structure of the generators is restricted in the fol- 
lowing way: For CSS codes, there are two types of gen- 
erators corresponding to Z- and X-type. We denote by 
ai a corner of the cube of generator type and by a'^ the 
body-opposite corner as depicted in Fig. [l] For non-CSS 
codes, the generators should satisfy ai — c^i- 

The number m of qubits per site should be bounded by 
the number of types of generators. For a local stabilizer 
code in any finite lattice, with open or periodic boundary 
conditions, there is a tradeoff in 3D 13^ 

kd = O(L^) 

between fc, the number of logical qubits, and d, the code 
distance, where L is the linear size of the lattice. If there 
are t < m types of generators, the number of independent 
stabilizers is at most tL^, and k > {m — t)L^. The code 
distance is then a constant independent of L. In order to 
achieve macroscopic code distance, it is mandatory that 
m < t. 

A string operator may wrap around a finite periodic 
lattice many times that it looks like a surface. But this 
is a property of boundary; as long as thermal stability is 
concerned, we ignore the boundary effects and consider 
stabilizer codes in the infinite lattice A = Z^. If a single 
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FIG. 1. Stabilizer Generators for non-CSS(top) and 
CSS(bottom) cubic codes. Throughout the paper we fix the 
coordinate system as shown. 



site operator i?, i.e., two qubit operator, is logical, we 
want it to be an element of the stabilizer group S. Since 
the stabilizer group does not explicitly include a single 
site operator, it is not easy to formulate the condition 
E £ S. For simplicity, we require that E is the identity 
up to phase. For a single site operator E, we denote by 
E[v]p the Pauli operator repeated along the line parallel 



to V passing p, i.e., 
E[v]p 



E®Ei 



(1) 



whose support, the set of sites on which a Pauli operator 
acts nontrivially, is the line, 

supp(£;[v]p) = {p + TO e A I n e Z}. 

We say E[v\ has period one if ||w||oo — 1- (||(ai^7c)||oo = 
max{|a|, |6|, |c|}.) We demand that any logical operator 
of period one be the identity up to phase. This condition 
is not sufficient for the code to be free of string logical 
operators, but is necessary. We will see that a nontrivial 
logical operator of period one is a string logical operator 
in our formal definition of strings in Sec. |V] 

Imposing the constraints above may result in a trivial 
code in a finite lattice for which the number of encoded 
qubits is zero (fc = 0). To avoid such a case, we restrict 
the generators such that the product of all corner oper- 
ators to be the identity operator up to phase. This con- 
dition is automatically satisfied by non-CSS codes under 
consideration. For CSS codes, this becomes a nontrivial 
algebraic constraint on the corner operators. 

There are equivalence relations on the set of codes. If 
two stabilizer codes are related by a symmetry transfor- 
mation of the unit cube, they are essentially the same. If 
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TABLE I. Complete list of cubic codes. The corners of the 
unit cube are labeled as in Fig[T] The second generator of 
non-CSS Code is given by the spatial inversion about body- 
center. The rest are all CSS codes, for which X-type generator 
is uniquely determined by eq. Q. The codes marked with ★ 
do not have string logical operators, while those with f do. 
See Theorem [2] and Appendix [B] 

one can be transformed into the other by a basis change 
on each site, we also regard them as the same codes. Re- 
naming of stabilizer generators obviously gives equivalent 
codes. Up to these equivalences, we report that there are 
1 non-CSS and 17 CSS cubic codes listed in Table [J The 
conditions of the cubic codes are summarized below. 

Condition 1: There are one or two qubits per site in 
the infinite simple cubic lattice Z'^. 

Condition 2: The stabilizer group S is translation- 
invariant and is generated by two types of opera- 
tors acting on eight corners of a cube. For non-CSS 
code, the two are related by spatial inversion, i.e., 
a'l — a2- See Fig. [1] The product of all corner 
operators of a CSS code is the identity. 

Condition 3: If e is a single site operator, then 
E is the identity up to phase. 

Condition 4: If ^ G has period one in Zoo-metric, 
i.e., supp(Z) is along one of 3 coordinate axes, 6 
face-diagonals, or 4 body-diagonals, then I is the 
identity up to phase. 

In the next section, we will study the conditions system- 
atically. 



III. COMMUTATION RELATIONS OF CORNER 
OPERATORS 

Given a set {gi, . . . , (7„} C Vm of n Pauli operators 
acting on m qubits, we can express their commutation 
relations in an n x n skew-symmetric (and, at the same 
time, symmetric) matrix u over the binary field. 

If we express g^'s in the columns of a 2m x n matrix P, 
then obviously 

P^AP = LO. 

Since lo is skew-symmetric, r = rank(a;) is even. 
Note that the rank of A is 2m. Since r < 
min{rank(P), rank(A)}, we see that r/2 is the minimum 
possible number of qubits on which Pauli operators of P 
act. Conversely, 

Lemma 1. Given a commutation relation uj of Pauli 
operators, all realizations Pofuj using minimum number 
of qubits are equivalent up to symplectic transformations. 

Proof. The rank of P is at least r = 2m. Being of full 
rank, P has linearly independent rows and we can add 
extra n — 2to rows to P so that the extension P'^ of P 
is invertible. Let Pi,P2 be two solutions realizing a;, 
and Pf , P| be their arbitrary invertible extensions re- 
spectively. We have 

(Pf )^A"P{= = w = {PifX'^P^ 

where 

Therefore, S"^ — Pf (PI is a symplectic transformation 
preserving A'^. The most general form of a transformation 
preserving A'^ is 

where S is such that A ~ S'^XS. Immediately, Pf = 
S^Pi, or Pi =SP2. □ 

We will translate all the requirements for the cubic 
codes into conditions on commutation relation matrix uj 
of corner operators. First, lo must represent a stabilizer 
code. The generators at different locations will commute 
if the components of uj satisfies a certain linear equation. 
We must consider all the cases when two generators meet 
with each other at a site, at an edge, at a face, and when 
they overlap completely. Note that this classification is 
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based on Condition 2. For non-CSS codes, the equations 
are 



uj{C,C') = 0, 



io{B,B') = 0, 
uj{D,D') = 



for the generator meeting at a site, 



Lo{A,C')+uj{C,A') = 0, 
uj{A,B')+uj{B,A') = 0, 
w(C,B)+w(B',C") = 0, 

for them meeting at an edge, 



uj{B,D')+uj{D,B') = 0, 
ujiC,D')+u;{D,C') = 0, 
uj{A,D) +uj{D',A') = 



uj{A, B) + cj(C, D) + A') + uj{D', C) = 0, 
uj{A, C) + uj{B, D) + uj{D\ B') + lo{C' , A') = 0, 
D') + C") + w(C, B') + wp, A') = 

for them meeting at a face. When generators meet at a 
cube, they automatically commute. Above 13 equations 
are independent since each term, e.g. u}{A,A'), appears 
only in one equation. For CSS codes, we only need to con- 
sider commutation relation between X-type and Z-type. 
There are 8 equations for the case when two generators 
meet with each other at a site, 12 equations when they 
meet at an edge, 6 equations when they meet at a face, 
and 1 equation when they overlap completely. One eas- 
ily checks that these 27 equations are independent since 
each term appears only once. 

Consider a non-CSS cubic code for which w is 8 x 8. 
Let m > 1 be the number of qubits per site, and P be the 
2m X 8 matrix whose columns consist of corner operators. 
A single site operator E is logical if and only if P^XE = 
0. Given E logical, we require = 0, or P^ X have the 
trivial kernel (Condition 3). In other words, rank(P) = 
2m. Therefore, we must have 



rank (io) — 2m. 



(2) 



Note that it implies m be the minimum possible realizing 
w. Conversely, if r = 2m, we need at least m qubits 
to realize oj, and rank(P) > 2m. Therefore, P^X has 
trivial kernel, and there is no nontrivial logical operator 
supported on a single site. 

Consider a logical operator E[y] of period one along 
y-axis. E ® E commute with A® B \i and only if 

X{A (E)B,E(E)E) ^ X{A, E) + A(B, E) = X{AB, E) = 0. 

See Fig. [I] Hence, E[y\ is logical if and only if 

X{AB,E)=Q X{CD,E)=0 



X{B'A',E) = 



XiD'C',E):^0 



These equations form a system of linear equations with 
unknown E, a 2m-component column vector. The coef- 
ficient matrix M is a 4 x 2m matrix. 

M = RP'^X 



where A is expressed in the first row of P^ , 
second, etc, and 



R = 



/^l 1 o\ 

1 1 

1 1 

yo 1 1/ 



B in the 



(3) 



Given ME — 0, we require E = (Condition 4). By 
eq.([2]), we see that P can be extended to be invertible, 
so that rank(i?P^A) = rank(i?(P^)'^A'=P'=) = rank(Pa;). 
The requirement becomes a simple formula: 



rank (Ruj) — 2m. 



(4) 



Conversely, if rank(Pa;) = 2m, then there is no logical 
operator of period one along y-axis. In an analogous 
manner, we consider all 13 logical operators along lines 
that are respectively parallel to 3 coordinate axes, to 6 
face-diagonals, and to 4 body diagonals. Note that the 
'derived' matrix Ru is calculated by adding rows of w cor- 
responding to corners that the logical operator of period 
one passes through. 

A CSS cubic code has 16 corner operators. The corners 
belonging to one of generators automatically commute 
with each other. Therefore, w has non-zero elements in 
off-diagonal blocks if we order the corner operators as 
{Az, Bz, . . . , Ax,Bx, ■ ■ ■}■ 



uj' 




Let us also order the basis of Pauli group "P,, 
Z-operators come first and 



such that 




The triviality of a single site operator is expressed as 
rank(a;) = 2m, or rank(w') = m. Consider X-type logical 
operator x[y\. It is logical if and only if 



X{AzBz,x) = Q 
X{B'zA'z,x)^0 



X{CzDz,x)^0 
X{D'^C'z,x)^Q 



which is equivalent to a matrix equation 



ppjx = 



(5) 



where x is an m-component vector, Az is expressed in 
the first row of Pj, etc, and R is given by eq.([3). Since 
u)' = {Pz)^Px bas rank m, Pz and Px are both of full 
rank m, and we can extend them to be P|, Px that are 
invertible. Since 



rank(i?Pj) = rank i?(P|)'^ 
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the matrix equation ([s]) is equivalent to 
m — rank(i?a;'). 



(6) 



As in the non-CSS case, there are 12 more equations 
ensuring the triviality of the X-type logical operator of 
period one. For Z-type logical operators, the equations 
are of form 



rank(i?a;'-' ) 



(7) 



We point out that it is a property of w whether or not 
the product of all corner operators yield the identity by 
Lemma[T] We have shown that the triviality of the single 
site operator (Condition 3) implies that any cubic code is 
a minimal realization of its commutation relation matrix 
w. Therefore, any two cubic codes with the same a; are 
related by a symplectic transformation, which is, in par- 
ticular, an invertible linear map. The product of all the 
corner operators of one code is zero (i.e., the identity), if 
and only if the product of all the corner operators of the 
other is zero. 

We thus completed the translation of the conditions 
for cubic codes into those for the commutation relation 
matrix of the corner operators. An advantage of this ap- 
proach is that we are classifying the codes up to symplec- 
tic transformation on sites. Moreover, the cases m = 1 
and m = 2 are treated simultaneously. 

There are 2"*^^^ w's of CSS codes that are consistent 
with the condition that the generators define a stabi- 
lizer codes, and 2""°"-°^^ cj's of non-CSS codes, where 
"CSS = 8 • 8 - 27 = 37 and nnon-css = 8*^2 - 13 = 15. 
After exhaustive search, we found no instance of a; sat- 
isfying the conditions when m = 1. Up to the symmetry 
group of the unit cube and renaming of generators for 
CSS codes (w' O w'^), we finally obtain Table [l] Here- 
after, we will call each code by Code 0, Code 1, etc., 
according to Table |lj 

The generators of CSS cubic codes show additional 
symmetries that we did not impose. Namely, Z-type 
generators and X-type generators are related by spatial 
inversion. Recall that ax, Q;^, a^, a'^ (a — A, B,C, D) 
denote the corner operators, each of which is a 2-qubit 
operator. Since they are purely Z- or X-type, we ex- 
press them by 2-component binary column vectors. For 
example, Az = ZI = (1 0)"^ of Code 1. We observe the 
following rule: 




ax 




(8) 



for all a = A, B,C, D. Because of this rule, there is a 
duality between X- and Z-type logical operators. That 
is, given the fixed origin of the lattice, for every X-type 
logical operator O, there exists a unique Z-type logical 
operator obtained by the spatial inversion about the ori- 
gin followed by the symplectic transformation on each 
site defined by eq.Q. Hence we wiU consider only X- 
type logical operators and Z-type stabilizer generators 
for CSS cubic codes. 



IV. MACROSCOPIC CODE DISTANCE 

In this section, we prove 

Theorem 1. Let d he the code distance of Code 0,1,2,3, 
or 4 defined on the periodic finite lattice . Then d > L. 

We introduce an important technique to deform a logical 
operator of cubic codes, which will prove the theorem. 
The technique depends on eq.(8|. We say a Pauli oper- 
ator is finite if its support is bounded (i.e., finite set). 
We will prove that any finite logical operator is a prod- 
uct of finitely many stabilizer generators. This implies 
Theorem [T] by the following lemma: 

Lemma 2. LetC{L) be a translation- invariant local sta- 
bilizer code of interaction range r > 1 (i.e., each genera- 
tor is contained in an r^ hypercuhe) defined on a lattice 
C^l)^ with periodic boundary conditions, where D is the 
dimension. Let d = d{L) be the code distance of C{L). If 
there exists Lq such that d{LQ) < Lg/lr — 1), then there 
exists a finite logical operator in the infinite lattice that 
is not a product of finitely many .stabilizer generators. 

We need a notion of connectedness: 

Definition 1. A set of sites {pi,p2, . . . ,Pn} is a path 
joining pi and p„ if for each pair {pi,pi^i) of consecutive 
sites there exists a stabilizer generator that acts nontriv- 
ially on the pair simultaneously, for i — 1, . . . ,n — I. A 
set M of sites is connected if every pair of sites in M are 
joined by a path in M. A connected Pauli operator is a 
Pauli operator with connected support. 

For example, {(0, 0, 0), (1, 0, 0)} is connected with respect 
to Code 0,1,2,3,4. {(0, 0, 0), (1, 1, 1)} is connected with 
respect to Code 2, but not connected with respect to 
Code 0,1,3,4. See Fig. [2j If a Pauh operator O is logical 
then any of its connected component is logical. If O is 
nontrivial, at least one of its connected components is 
nontrivial. 

Proof, (of Lemma [2]) Suppose d = d(io) < ^^/{r — 1) 
for some Lq. There exists a connected nontrivial logical 
operator O of support M, where the number of sites in 
M is d. Given a closed path {pi,p2, . ■ . ,Pn,Pn+i = Pi} C 
M, the union of the shortest line segments Ci : [0, 1] — 
connecting Ci(0) = pi and Ci(l) — Pi+i is a trivial 
homological cycle of the Z?-torus Z) {'Llo)'^ ■ If it is 
not the case, since any nontrivial homological cycle of 
has length Lq, we must have d{r — 1) > Lq. 
Consider a lifting of all closed paths of M into the uni- 
versal covering MP of via line segments cts. Since 
any closed path in M can always be express by a triv- 
ial homological cycle, the lifting is bounded. The corre- 
sponding lifting of O is not a product of finitely many 
stabilizer generators, since it was not trivial. □ 

The converse of Lemma [2] could be an interesting prob- 
lem, since it, if true, implies that the code distance of 
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FIG. 2. Stabilizer generators for non-CSS Code 0, and CSS 
Code 1, 2, 3, and 4. They all have code distance > L (Theo- 
rem [T|. The bottom four are free of string logical operators. 
See Section |V] 



translation-invariant local stabilizer code in periodic fi- 
nite lattice is either 0(1) or Q,{L). 

Consider a finite logical operator O. We will show 
that O is a finite product of stabilizer generators, i.e., a 
trivial logical operator. We may assume O is supported 
on a finite box _B C A, where 

B = z) \ xa < X < xi, yo <y <yi, zq < z < zi}. 

We first deal with Code 0. The stabilizer generators are 



depicted in Fig. [2j Consider the vertex v = (cci, yi, zq) of 
B that has largest x- and y-coordinate and the smallest 
z-coordinate. It must commute with ZX of Qo XY 
of g^. Since ZX ^ (1001) and XY = (0111) are inde- 
pendent, the commutation gives two constraints on v. A 
possible w is a linear combination of ZX and XY. (Re- 
call that Pauli group is abelianized to be a vector space.) 
li V = II, then we can shrink B, the support of O. If 

V = ZX, we can multiply Qq inside B to make v = II, 
hence shrink B. li v = XY, then Qo inside B will make 

V = II. liv = Y Z, then the product QqQq inside B will 
make v — II. In short, we have deformed the support B 
of O such that B now consists of one less site. See the 
second figure of Fig. [Sj 

The process can be done arbitrarily many times as long 
as the deformed B can contain a unit cube so that the 
multiplication by Qq , Qq or both only affects the sites in 
B. Since we started with the finite box, we end up with 
a support consisted of three thin rectangles (the third of 
Fig. [3]). To be precise, a thin rectangle Rt perpendicular 
to i-axis means the set of sites 

Rx = {(a;o, y,z) \yo<y <yi, zq < z < z^}, 

Ry = {{x,yo,z) \ Xq < X < Xl, Zq < z < Zl}, 

Rz = {{x,y, 21) I xo < a; < Xl, yo < y < yi}. 



of Ry that is not 



Consider the vertex v' = {xi,yQ,zo 
contained in the other two thin rectangles. It must com- 
mute with ZX, XZ of Q and XY, ZY of . Therefore, 
v' = // (the fourth of Fig. [3|. Continuing, we deduce 
that whole rectangle Ry \ {R^ U i?^ ) is the identity. Note 
that this procedure was possible because we were able to 
find an edge that has "sufficiently independent" corner 
operators. We call an edge is good for erasing if the argu- 
ment above works. Similarly, one can show that the other 
two rectangles R^ and Rz are also the identity. (The edge 
corresponding to XX — ZX of Qq and ZZ — XY of Qq 
is good for erasing, etc.) Thus, we have shown that by 
multiplying appropriate stabilizer generators inside B, 
we get the identity operator. 

Secondly, let us show that Code 1 has macroscopic 
code distance. It suffices to consider X-type logical op- 
erators. Let i? be a finite box that supports an AT-type 
logical operator. Consider the vertex v = {xi,yi,zo) of 
B that has largest x- and y-coordinate and the smallest 
z-coordinate. v commutes with IZ of Qf , and hence is 
either II or XI. If v = II we shrink B by one site. If 
V = XI we multiply Qi inside B to erase v. We again 
end up with three thin rectangles. Consider the vertex 
v' = {xi,yQ, zq) on the rectangle Ry perpendicular to y- 
axis that is not contained in the other two rectangles. It 
commutes with IZ, ZI of Qf . Therefore v' = II. Con- 
tinuing, we erase Ry. Similarly, one can erase the other 
two rectangles. Note that for CSS cubic codes, an edge 
is good for erasing if the corner operators of at the 
ends of the edge are independent. 

This strategy is good enough to show that the code dis- 
tance is macroscopic for Code 2, 3, and 4. We summarize 
the erasing procedure in Fig. [3j 
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\ Code 0, 1, 2, 4 




\ Code 3 

FIG. 3. Proof of macroscopic code distance. Deformation of 
a finite logical operator is depicted for each Code. For Code 
1,2,3,4, the logical operator of X-type is considered. 

V. LOGICAL STRING SEGMENTS 

A string logical operator might be regarded as a logical 
operator whose support is one dimensional. Indeed, the 
logical operators of some codes have definite topological 
structure. For the 2D Ising model or the toric code [U 
III HI] , the syndrome corresponding to a single site error 
has particular shape, by which we endow the distribution 
of Pauli operators with topological meaning. Concretely, 
X-error on 2D Ising model can be represented by a square 
in the dual lattice, and Z-error on 2D toric code by a link 
in the real lattice. 

However, the topological meaning of an operator may 
not always be well-defined. The most important prop- 
erty of string logical operators would be that it can be 
extended in the infinite lattice to a arbitrarily long string 
with constant width. It is in fact the property that is 
used in no-go theorems on quantum memory based on 
stabilizer codes in 2D [3H5]. We capture this property of 
string logical operator in the definition as follows. 

Definition 2. Let ili,f72 C A be congruent cubes con- 
sisting of sites, and O be a finite Pauli operator. A 
triple C = (O, ill, 1)2) is a logical string segment if ev- 
ery stabilizer generator that acts trivially (by identity) 
on both ill and ft2, commutes with O. We call r2i,2 the 
anchors of The directional vector of ^ is the relative 
position of ili to fl2- The length of C is the Zi-length of 
the directional vector, and the width is w. 

Since fli and Q2 are congruent, the directional vector 
of a logical string segment is well-defined. If the direc- 
tional vector is (a, b, c), its /i-length is laj-t-lfol-f |c|. As an 
example, a string operator of 2D toric code that creates a 
pair of vortex excitations is a logical string segment with 
anchors being the plaquettes carrying the vortex. Note 
that O may not commute with all stabilizers that acts 
trivially on the anchors. For O with supp(O) contained 
in the anchors is obviously a logical string segment. We 



need to exclude such a trivial case. 

Definition 3. A logical string segment ( = {0,rti,rt2) 
is connected if there exist two sites pi S ^i,P2 G ^^2 that 
can be joined by a path in supp(O) U {pi,P2}, where 
supp(O) is the set of sites on which O acts nontrivially. 
Two logical string segments {0,fti,fl2),{0',ili,il2) are 
equivalent if O' is obtained from O by multiplying finitely 
many stabilizer generators. C is nontrivial if every equiv- 
alent logical string segment is connected. 

For finite w, define ^(w) to be the maximum length 
of all nontrivial logical string segments of width w. 4> 
is a non-decreasing function on the set of positive inte- 
gers. Obvious from the definition is that > 3(w— 1) 
for any stabilizer code in 3D since any logical string seg- 
ment with overlapping anchors is always connected, and 
hence nontrivial. We allow <p{w) to assume infinite value. 
For example, </> becomes infinite at small values of w 
for 2D toric code [14], 3D toric code [7], and Chamon 
model [SI [T3] . String logical operator is defined as a logi- 
cal operator in the infinite lattice containing an arbitrar- 
ily long nontrivial logical string segment. 

Definition 4. A translation-invariant stabilizer code de- 
fined by a set of local stabilizer generators, is free of string 
logical operators if the maximum length (t>{w) of nontriv- 
ial logical string segment is finite for all finite w. 

2D Ising model is free of X-type string logical operators 
according to our definition. Consider an X-type logical 
string segment C — {0,fli,il2)- Being finite, O cannot 
be supported outside the anchors. Therefore, 4>{w) = 
2{w - 1) < 00. 

Theorem 2. Code 1, 2, 3, and 4 are free of string logical 
operators. 

We will argue that if the length of a logical string seg- 
ment is sufficiently larger than its width, then there exists 
an equivalent logical string segment that is disconnected. 
The rest of this section is the proof of Theorem [2] 

A. Reduction to flat segments 

Let C be a logical string segment of width w. Using the 
technique introduced in Section|IV| we will deform C such 
that it is a union of at most three logical string segments 
whose directional vectors are (a, 0, 0), (0, 6, 0), (0, 0, c) re- 
spectively; {a,b,c) is the directional vector of (. We as- 
sume C is of X-type. 

Consider C, of Code 1 such that a,b,c > 0. Since the 
support of C is finite, we can shrink it by multiplying 
stabilizer generators until C is contained in the smallest 
box B containing the anchors of C,. Note that this is 
possible due to a special property of Code 1. Namely, 
using the technique of Section [TV] one can shrink the size 
(> 1) of a finite logical operator by 1 from any direction, 
since there are two orthogonal edges that are good for 
erasing on each of six faces. 
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FIG. 4. Deformations of X-type logical string segments of 
Code 1 and 2. The small cubes and the filled dots are the 
anchors. The good edges for erasing are easily identified from 
the figures. 



Now the two anchors are located on (0, 0, 0) (after shift 
of the origin) and (a, 6, c), and the B is of size (a + w) x 
{b + w) X (c + w). Using the fact that ZI — IZ along x- 
axis of Qi (Fig. [2| has two independent corner operator, 
we shrink B as shown in the second figure in the first 
row of Fig. |4j Similarly, one shrinks B further using 
IZ— ZI (read rightward) along y-axis and ZI—IZ (read 
downward) along z-axis. The initial ^ has been deformed 
such that it is a union of three logical string segments 
that are parallel to the coordinate axes. One can easily 
extend the argument to the case where a, 5, c > 0. 

Observe that Code 1 has three-fold rotational sym- 
metry about (1, 1, l)-axis. Therefore, we only need to 
consider one more case where a,b > and c < 0. The 
strategy is the same as before. We finish the reduction of 
logical string segments of Code 1 by drawing the second 
row of Fig. |4] 

Code 2, 3, and 4 are treated similarly. For Code 2, the 
initial reduction of an arbitrary logical string segment to 
the smallest box that contains the two anchors is possible 
because there is a pair of orthogonal edges that are good 
for erasing on each of six faces of Qf . The subsequent 
reduction to flat segments is depicted in Fig. 4] Code 3 
and 4 needs more explanation. See Appendix X] 

Next, we show that each flat logical string segment is 
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FIG. 5. Deformation of flat logical string segement . The 
bottom two figures are implied by 'confusing constraints'. If 
C^y^ is long compared to its width, it is equivalent to a dis- 
connected one. 



equivalent to a disconnected one if it is long relative to 
its width. 



B. Confusing constraints 



Consider an JT-type logical string segment c}y'' — 
(O, ili,r22) of width w of Code 1 whose directional vec- 
tor is (0, Z,0), where I is the length of C^^', and O is 
supported on the smallest box B that contains f2i,r22- 
Place Qf such that it touches exactly two consecutive 
sites pi = (x, y, z), p2 = (a;, y -f 1, z) G B \ (VL^ U VL2) 
where x is the largest and z the smallest. Qf acts on 
V\ ~ Vi by IZ — ZI^ which gives a constraint on possible 
Pauli operators on pi — -pi since O is commuting with 
Qf . Explicitly, pi — p2 is a linear combination (with co- 
efficients in F2) of // - IX, XI - II, and IX - XI. We 
observe that given any possible operators on pi — p2 one 
can make them // — // by multiplying Q^'s inside B. It 
is an important and common property of a good edge for 
erasing (currently, it is IZ — ZI of Qf ), which is derived 
from eq .([8|) . It is a variant of the technique described in 
Section IrVl 

We shrink the support _B of O such that now B consists 
of the anchors plus two thin rectangles Rz , Rx where Ri is 
perpendicular to i-axis as shown in the top right figure of 
Fig. [5| Let us calculate what Pauli operators are possible 
along the edge e of R^ that is not contained in Ry. Any 
consecutive pair of two points on e away from the anchor 
commute with IZ — ZI and // — IZ of Qf . Hence, 
the pair must be a linear combination of XI — II and 
IX — XI. Since the base field is binary, there are only 
four combinations: 



II - II 
IX -XI 



XI -II 
XX -XI 



The edge e is a consistent sequence of such pairs. Go- 
ing to the right (increasing y-coordinate), // must be 
followed by //, XI by //, IX by XI which must be fol- 
lowed by //, and XX by XI which must be followed by 
//. That is, e is the identity except possibly 2 sites near 
the left anchor! We say such a constraint is confusing 
whose solution is eventually II. We can repeat the argu- 
ment on the next 'line' e' that has 1 smaller x-coordinate 
than e, to deduce that e' is the identity except possibly 
4 sites. Inductively, we conclude that Rz is the iden- 
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FIG. 6. Deformation of flat logical string segement Q . 

tity except possibly a small 'triangle' adjacent to the left 
anchor. See the bottom left figm^e of Fig. [5] 

We do the same calculation for Rx with constraints 
given by IZ — ZI and // — IZ from the bottom face of 
Qf . The algebra is the same. Summarizing, we have 
shown that if / > w + 2{2w — 1) = — 2, then C,^^ 
is equivalent to a disconnected one. Note that Code 1 
has three-fold rotational symmetry x — > y — )■ i — )■ i. 
Therefore, the maximum length 0i(w) of nontrivial log- 
ical string segment of Code 1 satisfies 

< \hw — 6 < oo, 

which completes the proof of Theorem [2] for Code 1. 

C. Inconsistent quasi-period 

Code 2 exhibits no three-fold symmetry, but instead 
two-fold symmetry (i O y) about the plane of normal 
vector (1,-1,0). Therefore, it suffices to consider two 
logical string segments C,y along y-axis and C,z along z- 
axis. 

Let (0,0, Z^) be the directional vector of Q ^ — 
(Oz, iliz, 5122) pertaining to Code 2, where Oz is sup- 
ported on the smallest box that contains the two anchors. 
Using the edge ZZ — Zl (read downward) of Qf , which 
is good for erasing, we further deform O similar to Fig. [5] 
The support B of is now the union of two rectangles 
Rx, Ry Q2 acts on the edge e of i?^: that is not contained 
in Ry by ZZ — ZI and IZ — ZI (read downward). This is 
a confusing constraint because any two consecutive sites 
on e must be one of 



II- II 
II -IX 



IX -XI 
IX - XX. 



Therefore, e is the identity possibly except two sites near 
the bottom anchor. The same inference is applicable to 
Ry due to the two-fold symmetry. The length of Q 
satisfies < w + 2w = Sw if ci^'' is nontrivial. 

Let (0,^y,0) be the directional vector of Q = 
(Oy, fiiy, ri2y). Using the edge ZZ — ZI of Qf along 
y-axis, which is good for erasing, we may assume that 
Oy is supported on the union of two rectangles Rz , Rx as 
in Fig. [6j The constraints on the outer edges ei of Rx 
and 62 of Rz (e2 has bigger a;-coordinate than ei) given 
by Q2 are not confusing. We need a different argument. 

Definition 5. A function / on the positive integers (?.e., 
a sequence) is eventually periodic with quasi-period t> 1 
and offset tiq > Q li f{n + 1) = f{n) for all n > uq. The 
period of / is the smallest quasi-period. 



If ti and t2 are quasi-periods of /, then f{n) = f{n + 
iti+jt2) for sufficiently large n where i, j are independent 
of n. Since there exist i,j such that iti + jt2 = t = 
gcd{ti,t2), t is also quasi-period. Therefore, 

Remark 1. The period divides any quasi-period. Note 
also that if two sequences /i and /2 have different periods 
ii,i2 and offsets ni,n2 respectively, there exists n' < 
max(ni,n2) -f lcm(ti,t2) such that fi{n') ^ f2{'n'). 

(2) 

Supposing Q is nontrivial and arbitrarily long, we 
will derive a contradiction: The operators on 62 has pe- 
riod 3 and quasi-period a power of 2. 

The constraints on a pair of sites of 62 are ZI — IZ 
and ZZ — ZI. Hence, the pair is one of 



II -II 
IX -XI 



XI- XX 
XX - IX. 



(9) 



The only possible infinite sequence on e2 is thus • • • — 
XI — XX — IX — XI — • • • , whose period is 3, or • • • — 
II — II — ■ ■ ■ . If 62 is the identity, then we consider next 
rows of Rz until we get a nontrivial row. If the entire Rz 
is the identity, then imposes constraints on lower edge 
of Rx, ZI — IZ and IZ — ZZ, whose solutions are given 
by eq.([9|. The period 3 is again revealed if nontrivial. 
Consider constraints on ei, or more generally, on four 
sites of Rx that form a square. 



-ij" 



IZ - ZZ 



ZI- ZI 



ZI - IZ 



ZZ - ZI 



Here we denoted each site operator as an element of a 
vector-valued matrix a. Note that, for example, XI and 
IZ is commuting because (10) (01)"^ = 0. Similarly, the 
fact that Rx is commuting with Qf can be expressed by 
a system of homogeneous equations: 

fli-lj-l-l + I „ , I ^-i-ij' 





or 




ai,j = 0, 



^(aj-ij'+i, ij-ij')- 



(10) 



In order to find possible operators on ei, we set ooj' — 
for aU /. Then, 



It is obvious that {aij}j is eventually periodic with pe- 
riod 1. For the rows with bigger i, the following holds: 
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Lemma 3. Suppose a set {ai_j}i_j of vectors over the 
binary field ¥2 satisfies 

Oij+i — Moij + 6(ai_ij+i, Oi-ij), 

where M"^ = M . If {ai_i^j}j>i is eventually periodic 
with quasi-period t and offset Uq, then so is {aij}j">i 
with quasi-period 2t and offset hq + t. 

Proof. Consider a sequence {cj}j>i given by 

Cj+i — Mcj + c', 

where c' is a constant 2-component vector over F2 inde- 
pendent of j. 

Cj+3 = Mcj+2 + c' = M{Mcj+i + c') + c' 

= Mcj+i + Mc' + c' = Mcj + Md + Afc' + c' 
— Mcj + c' = Cj_|_i 

Therefore, {cj}j>i is eventuahy periodic with quasi- 
period 2 and offset 1. Define 



where j > 1. Clearly, Ci_ij+i = aj_i^„„+i+jf = 
ai_i,„o+i = Ci_i,i, i.e., {ci_ij}j>i has period 1 with 
offset 0. Moreover, {cij}ij satisfies 

= Mci^j -\- c'i j, 

where c'^ j = d (oj-i ; 4- 1 -f- (j - l)t < / < no -h 1 jt) 
does not depend on j because {a^-i j}j is eventually peri- 
odic with quasi-period t and offset tiq. Therefore, {q 
is eventually periodic with quasi-period 2 and offset 1. 
Since „(,_|_;j_|_(j_i)4(l <h<t)is determined by Cij and 
{ai-i^no+h'+{j^i)t}i<h'<h, we get the claim. □ 

We have shown that each row of that is not con- 
tained in Rz is eventually periodic, and as we decrease by 
1 the z-coordinate of a row of R^ we get its quasi-period 
doubled. Since the initial row has period 1, the quasi- 
period of each row is some power of 2. We need to show 
the same thing continues to hold as we move toward 62. 
The constraints pertaining to the last edge of R^, i.e., 
the intersection of Rx and Rz , are 



IZ - ZZ 



II -II 



ZI - ZI IZ- ZZ. 

The recursive equation is then 

)• (11) 

One can repeat the proof of Lemma [3] except that now 
one has to show that a sequence defined by 

= Nf, + /', 



where iV^ = / and /' is a constant vector, is always 
eventually periodic with quasi-period some power of 2. 
This is easy: 

/,+2 = N{Nf, + /') + /' = /, + /" 

/,+4=/,+2 + = +/"+/" = /,. 

We have proved 

Lemma 4. Suppose a set {ajj}jj of vectors over the 
binary field ¥2 satisfies 



1 = Na.^j -(- 6(aj_ij+i,ai_i, 



where N"^ = /. //{ai_ij}j>i is eventually periodic with 
quasi-period t and offset uq, then so is {a-i,j}j>i with 
quasi-period At and offset ng . 

The exact recursive equation for the next row is slightly 
more complicated since there are stabilizer generators 
meeting three rows, {w — 1), w, (w -\- l)-th. However, we 
do not need detailed information how acts on (w — l)- 
and w-th rows to infer the quasi-period of (w-f l)-th. The 
constraints on i-th row {i > w) are as the following. 



ZI -IZ ZZ - ZI 

IZ - ZZ ZI - ZI. 

The recursive equation is then 



a 



1 0' 

1 1, 



V«i-2j + l, a,i-2.j , 



where i > w. We may regard that both (i — 2)-th and 
{i — l)-th rows have the same quasi-period and the same 
offset. Hence, we can apply Lemma |4] to the rows in Rz. 

We have reached 62 starting from ei. Since ei has 
period 1, and each following row has quasi-period some 
power of 2, 62 must have quasi-period 2^~'^ x4^ = 2^^~'^ 
with offset 2"^"^ — 1. The true period 3 of 62 must di- 
vide 2^™^^, which is a contradiction. By Remark [l] we 
conclude that the length of a nontrivial logical string seg- 
ment ^ is < w-f (2^-2 _ 1) + 3 . 23"'-2. Since the length 

of Q is < 3w, the maximum length 4>2(w) of nontrivial 
logical string segments of Code 2 satisfies 



4>2{w) < 5w + 2' 



< 00, 



which completes the proof of Theorem [2] for Code 2. 

The proof of Theorem [2] for Code 3 and 4 using similar 
technique is given in Appendix \K\ 
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Code fc(L) 

L + 3 ■ 2" (92 + 2q7 + 8(79 + 48(763 + 64g65 + ISqm) 

1 2 [1 - 2g2 + 2''+^ (ga + 12gi5 + GOgea)] 

2 2''+i (1 + 6(77 + 6(721 + 30g3i + 60(763 + 126(7127) 
g 2'-+i(l + 8(715 + 6921+ 40^31 +42q63 

+ 16^85 + 112gi27) 

4 I 2'-+' (1 + 2(73 +8gi5 + 40(731 + 48(763 + 112gi27) 

TABLE II. Exact empirical formulae for the number k of en- 
coded qubits in periodic finite lattice as a function of L 
(2 < L < 200). Here, q„ = qn{L) is the divisibility function, 
and r — r{L) is the largest integer such that 2^ divides L. 




FIG. 7. Tunnel used to derive the upper bound on k, the 
number of encoded qubits. The stabilizer generators outside 
the tunnel are independent. 



Code 


Lower bound 


Upper bound 





L + 6q2 


12L- 12 (4L if TfL) 


1,2,3,4 


2 


4L 



TABLE III. Lower and Upper bound on the number of en- 
coded qubits in periodic finite lattice Zf,. Here, g„ = qn{L) 
is the divisibility function. 



VI. NUMBER OF ENCODED QUBITS AND 
LOGICAL OPERATORS IN FINITE PERIODIC 
LATTICES 

A. Number of encoded qubits 

The number of enco(ied qubits of a stabilizer code in 
any finite periodic lattice will be obtained, once we know 
all the algebraic relations of stabilizer generators in the 
infinite lattice. A nontrivial example of this approach 
is given in [15] for Chamon model. It has a nice prop- 
erty that the product of stabilizer generators becomes 
the identity only when they form body-diagonal surfaces 
in the infinite lattice. Since there are 4 body-diagonals, 
the number of encoded qubits of Chamon model in the 
periodic lattice Z2p^ ^^2py ^'^2p, is fc = 4:gcd{px,Py,Pz)- 

Our cubic codes exhibit even more peculiar dependen- 
cies of k on the linear size of the periodic lattice Z| . We 
found empirical formulae for k = k{L) of Code 0, 1, 2, 3, 
4, which are exact if 2 < L < 200. For ease of notation, 
we define the divisibility function Qn on positive integers 
for each positive integer n by 



QniL) 



1 if n divides L, 
otherwise. 



The formulae for k are given in Table |TTj 

Although we only know empirical formulae of k for 
small values of L, we can prove lower and upper bounds 
on k. Since there are equal number of stabilizer gener- 
ators and physical qubits, k is equal to the number of 
independent algebraic relations of generators. It is obvi- 
ous from the definition of cubic codes that k > 2 since 
the product of all generators is the identity in any peri- 
odic lattice. In addition, k is always an even number for 



CSS cubic codes because of the duality between X- and 
Z-type generators. 

The non-CSS Code is more complicated than the 
CSS codes. We prove the bounds on k for Code in 
Appendix [c| Consider Code 1 in We find an inde- 
pendent set of generators of Z-type to derive the upper 
bound. Let S be the set of Z-type generators lying out- 
side a straight tunnel T of length L parallel to z-axis, 
whose cross-section is an 1 x 2 rectangle enclosed by 6 
sites. There are — 2L generators in S. We claim that 
S is an independent set of generators, and hence k < 4L. 

Suppose a linear combination O of generators in S is 
the identity operator. We show O is the zero combina- 
tion. Choose the origin of the coordinate system such 
that the sites enclosing the cross-sectional rectangle of 
T are described by a; = ±1, y = 0,1. Let Z^^^ be the 
straight hue given by a; = 0,y = 1. See Fig. [Tj Since 
O is the identity, in particular, Z^^^^ is acted on by the 
identity. Every unit edge in l^^'^ connecting (0, 1, i + 1) 
and (0, l,i) is one of (read downward) 



II -II 
IZ -II 



ZI-IZ 
ZZ -IZ 



(12) 



which is canceled by the neighboring unit edges. If e\ — 
II - II, then €2 = II - II. If ei = ZI - IZ, 63 ^ 
IZ - II. If ei ^IZ- II, €2 = II- II. If ei^ ZZ- IZ, 
62 ~ IZ — II. We see that {ei}i>i is eventually II — II. 
Since Z'^^ is periodic, II — II = cl+i = ei. We conclude 
that the coefficients of generators in S around Z*^^^ are 
all zero. Now O is a linear combination of generators 
lying outside the enlarged tunnel whose cross-sectional 
rectangle is described bya; = ±l, ?; = 0,2. 

We can repeat the argument inductively on the lines 
parallel to z-axis, each of which is given by 



{iO,y',z)eA\0<z<L}, 



where y' = 2, . . . , L — 1, until the tunnel becomes a slab 
of width 2. (The slab is in between two planes x — ±1.) 
Consider a straight line Z'^^-' parallel to z-axis given by 
y ~ 0,x — 1. Since possible operators on each unit edge 
of r'^-' are again given by eq.(12|, we argue similarly to 



conclude all coefficients in O are zero. This completes 
the proof of the upper bound on k of Code 1. 
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For Code 2, we use the same initial tunnel T. The 
initial is chosen to be the line parallel to z-axis given 
by X = 0, y = 0. For Code 3, we choose l^^'^ — Every 
unit edge of l^^^^ is one of 



II - II 

IZ - zz 



IZ - II 

II - zz 



which must be canceled by neighboring unit edges. The 
only choice is // — //. Due to two-fold symmetry of 
Code 2 and 3, the rest of calculation is easy and proves 
the upper bound. For Code 4, we choose the same initial 
tunnel T and the initial line Z*-^' = Arguing as above, 
one enlarges the tunnel until it becomes a slab of width 
2. Now consider any horizontal line V {z = z') in a; = 1 
plane. We see that every edge is one of 



II -II 
II -IZ 



ZI -IZ 
ZI -II 



which must be canceled by neighboring edges. The only 
choice is // — //. This completes the proof for the upper 
bound on k of Code 4. 



B. Plane logical operators 

Theorem [2] says, in particular, that Code 1, 2, 3, 4 
do not have any nontrivial logical operators that is sup- 
ported on a thin strand {l,...,w} x {l,...,w} x 
where w — 0{\). Bravyi and Terhal [3] showed that 
any local stabilizer code must have nontrivial logical op- 
erator supported on a thin slab {1, . . . , w} x Z^, where w 
is the interaction range. Indeed, cubic codes have logical 
operators on (w = l)-slab, i.e., plane logical operators. 
We consider the simplest plane logical operators that are 
repetition of a single site operator (two qubit operator). 
For a single site (two qubit) operator we define 



E 



(a,b,c)-pla 



to be the tensor product of E over the plane orthogonal 
to (a, 6, c). cr|; is logical if and only if 

\{E, Az) + \{E, Bz) + KE, Cz) + \{E, Dz) = 0, 
\{E, A'z) + \{E, B'z) + \{E, C'z) + X{E, D'^) = 0. 

which is equivalent to \{E, AzBzCzDz) — and 
\{E,A'^B'2C'zD'^) = 0. (The two arc in fact equiva- 
lent because one of our conditions defining cubic codes 
requires that the product of all eight corner operators be 
//.) We see that Uj^"' is a logical operator of Code 1, 



[1,-1,0] 

zz 



'IX 

Moreover, if the linear lattice size L 



and so is a 

is odd, they anti-commute, and hence are both nontriv- 
ial. In this way, one can easily find logical operators of 
form cr^'"^'. We make it clear that the plain logical op- 
erators found in this way do not generate all the logical 
operators. Empirically, for some special lattice size, e.g. 
L = 8, 15, 63, there are many more logical operators that 
cannot be described by the plane logical operators. 



{XI) 



{IX) 



FIG. 8. Syndrome cubes of Code 1 caused by a Pauli operator 
of weight one. Filled dots indicate excited Qf . See also Fig. [2] 



VII. THERMAL STABILITY 

In this section, we discuss the energetics of implement- 
ing nontrivial logical operators of Code 1. A natural 
choice of Hamiltonian for a local stabilizer code is the 
sum of local generators: 



(13) 



peA 



whose ground space is the code space. We consider ad- 
verse logical operations on the code space by thermal 
environment. Quantum tunneling between two different 
ground states is exponentially suppressed in the system 
size since eq.(13) is gapped and the code distance is at 
least L. We may model the thermal noise as a Marko- 
vian chain of actions to the system by Pauli operators of 
weight one. At each step, a Pauli operator by the envi- 
ronment will excite some terms in the Hamiltonian. Such 
excitation is completely determined by the accumulated 
syndrome, the data that describes which stabilizers are 
unsatisfied. Formally, a syndrome can be viewed as a 
Z2-valued function on the stabilizer group. 

We claim that the energy of any partial implementa- 
tion of cr^j^^ is proportional to its boundary length. A 
partial implementation of a logical operator is a Pauli 
operator that is a restriction of the logical operator on 
a subset of sites. It has excitations, if any, only along 
the boundary. Note that a syndrome caused by a Pauli 
operator is the sum of syndromes caused by Pauli oper- 
ators of weight one, each of which is expressed as a cube 
in the dual lattice as shown in Fig. [8j The syndrome 
corresponding to the partial implementation of cr^^*^' is 
expressed as a stack of /X-cubes. An outer boundary 
point of the stack is either a vertex of one syndrome 
cube, a point where two cubes meet, or a point where 
three cubes meet. In any case, one can easily verify that 
there exists a filled dot i.e., an excitation within distance 
one from the boundary point. Therefore, the number of 
excitations (the energy) of any partial implementation of 

"'/x^' least the length of the outer boundary. One 
checks that the same holds for plane logical operators 

[110] , [1,-1-0] 

CT^/and(T^^ ■ 

This property that energy is proportional to the 
boundary length of partial logical operator mimics that 



13 



of 4D toric code. Still, it is needed to show that any 
partial implementation of any nontrivial logical operator 
confronts high energy barrier. But we do not have a com- 
plete description of all equivalent variants of, say, crf^"^ 
of Code 1. Moreover, we do not have a list of logical 
operators for arbitrary system size. 



We finally note that the system governed by eq.(13) 



is a quantum glass any isolated excitation (defect) 
cannot propagate easily to meet its partner to annihilate 
if the temperature is low. The situation is more stringent 
than Chamon model where one can move even number 
of excitations by string operators. In Code 1, there is 
no local Pauli operator that can move a localized, but 
not locally created, set of excitations to nearby location 
congruent ly. 



VIII. DISCUSSION 



We have defined logical string segments for local stabi- 
lizer codes, whose length determines the existence of non- 
trivial string logical operator. We classified translation- 
invariant CSS stabilizer codes in three dimensions with 
two stabilizer generators, and also found the unique non- 
CSS stabilizer code with the special symmetry. We 
showed that some of the CSS cubic codes do not have 
string logical operators. The codes without string logical 
operator exhibit peculiar dependence of the number of 
logical qubits on the system size. This is because there 
are complicated algebraic relations of stabilizer genera- 
tors, for which we do not have full description. 

The thermodynamic stability of encoded information 
in the code is still an open problem. It has been proved 
that in order for a quantum memory based on stabilizer 
codes to have a thermal stability, it suffices to have a 
good error correcting procedure by which every low en- 
ergy syndromes are good [16^ . Returning to Code 1 with 
the system size such that fc = 2, there is only one al- 
gebraic relation of stabilizer generators: The product of 
all generators is equal to the identity. Therefore, there 
is the unique constraint (conserved charge) on the space 
of syndromes: The number of excitations must always 
be even. In other words, any state with two excitations 
arbitrarily far apart is allowed. Since there is no string 
operator, a Pauli operator that causes such excitations 
must be of complex shape. A good error correcting pro- 
cedure should answer how such a syndrome is formed. 
On the other hand, if every three dimensional local sta- 
bilizer code is not stable, our codes illustrate that one 
should not attempt to prove it by showing the existence 
of string logical operator. 

It has been recently proved that any Pauli walk that 
results in a nontrivial logical operator of Code 1 must 
experience energy barrier of height ri(logL) |17j . 
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Appendix A: Proof of Theorem [2] for Code 3, 4 

For any logical string segment C, we first deform it such 
that it is supported on the smallest box containing two 
anchors. This is easily done if there are two orthogo- 
nal edges that are good for erasing on each face of 4. 
However, Code 3 and 4 have faces that do not have such 
two edges. Let II of Qf be at (0, 0, 0). The (x = 0)-face 
of Qf has only one edge IZ — ZZ that is good for eras- 
ing. Using this, we can shrink the x-size of the support 
of C from positive x-axis, except a y-plane Ry of thick- 
ness one. The constraints on the outer edge e of Ry are 
IZ — II, IZ — ZZ. Every unit edge in e is thus either 
// - //, XI - II, II - XX, XI - XX. Since e is finite, 
every site of e is //. The [y = 0)-face of Qf is sym- 
metric to (x = 0)-face. Code 4 has two faces that do not 
have two orthogonal good edges for erasing: {x = 0)- and 
[z = 0)-face of Qf. One can repeat the same argument. 

Fig. [9] shows that every logical string segment of Code 
3 and 4 is equivalent to a union of logical string seg- 
ments whose directional vectors are parallel to the co- 
ordinate axes. Code 3 has two-fold symmetry i; o y, 
and therefore it is enough to consider Q along z-axis 



and Cy'^' along y-axis. C 



(3) 



IZ - ZZ {read downward) of Qf 



can be further deformed by 
See Fig. 



10 



We con- 



sider the constraints IZ — II, IZ — ZZ given by Q3 on 
the external edge. They are confusing because we have 
to obtain a consistent solution on the edge using 



II -II 
II -XX 



XI -II 
XI - XX. 



Therefore, a nontrivial Q has length < w + {2w — 1). 

Cy'^' is deformed using ZZ — Z/(read rightward) of Qf 
as depicted in Fig. [TO] One edge is confused, and we are 
left with a rectangle normal to z-axis. ei is constrained 
by IZ — IZ and // — ZZ . A possible neighboring pair 
is one of // - //, XI - II, IX - XX, and XX - XX. 
Thus, any nontrivial ei is eventually • • • — XX — ■ ■ ■ of 
period 1. 62 is constrained by IZ — ZZ and ZZ — ZI. 
A possible neighboring pari is one of // — II, XI — XX, 

IX - XI, and XX - IX. A nontrivial 62 is XI - 

XX — IX — XI — ■ ■ ■ of period 3. In between ei and 62, 
we have 
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T 



















Code 3 




\ Code 4 

FIG. 9. Deformations of logical string segments of Code 3, 4. 
Anchors are marked by the filled dots. 



confused edge 



FIG. 10. Deformation of flat logical string segments Q and 

^(3) 

■ 



(read downward) are confusing constraints. must be 
short because 

{ZI - II, ZI - IZ} 
{ZI - II, IZ - ZI} 

(read as decreasing x-coordinate) are confusing con- 
straints, (^y'^^ is deformed to have eventually periodic two 
edges ei, 62. On ei, the constraint is {// — IZ, ZI — ZI} 
whose solution is a linear combination of IX — II and 
XI — XI. Hence, ei is of period 1. On 62, the constraint 
is {IZ — ZZ, ZI — IZ} whose solution is a linear combi- 
nation of IX — XI and XI — XX. Hence, 62 is of period 
3. In between ei and 62, we have: 



H-l,j 



ZI -IZ 



IZ -ZZ 



O-i-lJ Oi-lj + l 

II - ZZ IZ- IZ 

ZZ - ZI IZ - ZZ. 

The recursive equation is then 



// -IZ ZI - ZI 

The recursive equation is 

flij+i = ^ flij + 6(a,_ij, ai_ij+i), 

to which we apply Lemma [3] We completed the proof of 
Theorem [2] for Code 4. 



Appendix B: Cubic codes with string logical 
operators 



Since AI ^ Jj ^ quasi-period of z-th row 

is twice as that of (i — l)-th by Lemma |3] Therefore, 
a nontrivial (y^'' has finite length. See Remark [l] This 
completes the proof of Theorem [2] for Code 3. 

Code 4 has no symmetry, and we need to check 

C^^"* , Cy^^'X^^ ■ The calculation is straightforward foUow- 
Ci^'* must be short if nontrivial, because 



ing Fig. 



11 



{ZI - II,IZ - ZI} 
{IZ ~ ZI, ZZ - ZI} 



For a local stabilizer code defined on a periodic finite 
lattice Z|^, suppose there exists a nontrivial logical opera- 
tor O supported on a strand {1, . . . , w} x {1, . . . , w} xZ^. 
Ifu'-|-2(r— 1) < L where r is the interaction range of 
the stabilizer generators, then we may lift O to a periodic 
logical operator O' of the code on the infinite lattice sup- 
ported on {1, . . . , w} X {1, . . . , w} X Z. Suppose further 
that any finite logical operator is trivial in the infinite 
lattice. If a contiguous part C of O' is trivial as a logi- 
cal string segment, then every congruent part of C is also 
trivial; O' is a disconnected union of finite logical op- 
erators. Hence O is trivial, a contradiction. Therefore, 
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FIG. 11. Deformation of flat logical string segments ^i*', (^'^ 



and ^i;"*' . 



Code 


string logical operator 


complement (L — 5) 


11 


■^■2[2](000)-^-^[5](100) 


^-^[y]{ooo)-^^[j/](ioo) 


12 


-^^[z](000)^-'"[«](010) 


-'^-^'[^](000)-'f^[*](010) 


13 


Z Z [z] (^OQO) I Z[z] (^oio) 


4°/' (000) 


14 


IX[z](^ooo)XI[z](oio) 


4°f'(ooo) 


15 


Z I [y] (000) ZZ[y] (^100) 


-^^[z](000)^-^'[2](100) 


16 


ZZ[101](ooo)IZ[Wl](,oo) 


IX[110](ooo)^/[110]{ioo) 


17 


ZZ[x](ooo)IZ[x](ooi) 


4T'(000) 



TABLE IV. Nontrivial string logical operators. -Eftilp denotes 
the Pauli operator such that E is repeated along v passing 
through p, and (p) denotes the Pauli operator such that E 
is repeated on the plane containing p perpendicular to v. 



any contiguous part of O' is a nontrivial logical string 
segment. The maximum length of the nontrivial logical 
string segment is infinite. 

In Table |IV| we list nontrivial logical operators sup- 
ported on strands for the cubic codes marked with f 
in Table |l] Recall that E[v]p denotes the Pauli operator 
■ ■ ■ ^ E (S) E ^ ■ ■ ■ on the line along the vector v passing 
through p. For example, ZI[0, 0, 1] (3,0,0) ^^[0, 0, 1] (2,0,0) 
represents (when L = 5) 



// // ZI ZZ II 

II II ZI ZZ II 

II II ZI ZZ II 

II II ZI ZZ II 

II II ZI ZZ II. 



Appendix C: Code 



1. Three- fold symmetry 



We remark that there is a three-fold symmetry for the 
generators of Code 0. (Fig. [2| If we rotate Qq by 120° 
about (1, 1, l)-axis, and then apply the transformation 



IX -^IZ — we see that Qo is invariant. This 

is in fact expected from the commutation relation w of 
corner operators of Qo; "we calculated Qo from lo. The 
trasformation is symplectic. If we order the basis for 
V2 as {XX., XZ^ ZX, ZZ}, the transformation is S'i?i2o° , 
where -Ri2o° is the rotation about (1, 1, l)-axis and 

/l o\ 
10 
11 

yo 1 0/ 

One directly checks that 5 preserves the symplectic form 
of the abelianized 7^2 

/O 1 o\ 
10 
1 
yo 1 0/ 

i.e., XS = X. This is an example of Lemma[l] A direct 
consequence of the three-fold symmetry is that logical 
operators always appear as a triple. 



2. Bounds on the number of encoded qubits 

We show fc, the number of encoded qubits of Code 
defined on the periodic finite lattice Z^, satisfies k > L. 
To this end, we present an algebraic relation of stabilizer 
generators in the infinite lattice that can be embedded 
into a periodic lattice of arbitrary linear size. An al- 
gebraic relation of stabilizer generators is a distribution 
of generators whose product is the identity. (Formally, 
the space of algebraic relations is the kernel of the linear 
map {f : A* G} ^ {(7 : A ^ V2} where A* is the 
dual lattice to A = Z'^, G is the abelian group of labels 
of stabilizer generators, and V2 is the abelianized Pauli 
group on a single site. G is isomorphic to Z2 x Z2 for 
our cubic codes because there are two types of genera- 
tors. This map is meaningful for translation-invariant 
stabilizer codes.) Note that the set of the locations of 
generators in an algebraic relation need not be finite; the 
product is well-defined if the number of generators acting 
on each site is finite, i.e., locally finite. For some relations 
the set of the locations of generators form a sublattice, 
which we call relation lattice, and can be described by 
a unit cell and basis vectors. Indeed, the relations that 
gives the lower bound k > L forms sublattices. 

The relation we consider first has its own lattice struc- 
ture i?i with a basis {(1, 0,-1), (0, 1, —1)}, the unit cell 
of which consists of Qq and Qq such that Qq is at (0, 0, 1) 
relative to Qq. In other words, Fig.[T2]is repeated accord- 
ing to (1, 0, — 1), (0, 1, — 1). Since the two basis vectors 
have period L in the finite lattice of linear size L, one 
can embed this relation into any finite periodic lattice, 
and there are L linearly independent such embeddings 
via translations along (1,0,0). Therefore, we have at 



16 



XY 



ZY 




combination of 



ZY 



XY 



FIG. 12. Unit cell of a relation lattice of Code 0. The relation 
lattice has basis {(— 1, 0, 1), (0, — 1, 1)}. The boxes displays 
the configuration of Qo's in the relation. Qq lies on the top 
of Qo. The line designated by the arrow is acted on trivially 
by four unit cells around it, as one can directly check using 
the Pauli operator diagrams. 



least L independent algebraic relations of stabilizer gen- 
erators for a L X L X L lattice. This proves that k > L. 
Note that this relation lattice is invariant under the three- 
fold symmetry of Code 0. If we consider the embeddings 
into finite periodic lattice with arbitrary three linear sizes 
Zl^ X Zl^ X Zl^, we see that k > gcd{L.j;, Ly, L^). 

Another relation lattice i?2 has basis 
{(2, 0,0), (0,2,0), (0,0, 2)} with unit cell such that 
Qq is at (1,1,1) relative to Qo- R2 is embedded 
into any finite lattice of even linear size. There are 8 
independent relations whose underlying lattice structure 
is i?2- However, i?i and R2 are not all independent. It 
is not hard to see R2 gives 6 more algebraic relations. 
This proves the lower bound k > L + 6(72 ■ 

It should be pointed out that — / ^ S is not auto- 
matically guaranteed since Code is non-CSS. However, 
the problem is resolved if one choose appropriate sign for 
each dependent generator. One can check directly that 
relations Ri and R2 does not generate — /. 

In order to prove an upper bound on fc, suppose 7 does 
not divide L. Then we can show that k < AL. The proof 
is very similar to that of Code 1. Let S be the set of 
generators lying outside a straight tunnel T of length L 
parallel to y-axis, whose cross-section is 1 x 2 rectangle 
enclosed by 6 sites, i.e., the tunnel is 1 x L x 2 {Ax x 
Ay X Az). There are 2L^ — AL generators in S. We show 
S is an independent set of generators. 

Suppose a linear combination O of generators in S is 
the identity operator. We show O is the zero combina- 
tion. Choose the origin of the coordinate system such 
that the sites enclosing the cross-sectional rectangle are 
described by x = 0, 1, z = —1, 0, 1. Let I be the set of 
sites on the straight line given by a: = z = 0. Since O 
is the identity, I is acted on by the identity. Every unit 
edge Bi in I connecting (0, i, 0) and (0, i -I- 1, 0) is a hnear 



XX - ZI 
ZY -XY 



ZZ -II 

XZ-ZX (CI) 



which is canceled by the neighboring edges. Since the 
left-hand side operators of eq.(Cl| are independent, we 
can unambiguously determine e^+i given e^, e.g., if ei — 
XX - ZI, then e2 = ZI - XZ. The right-hand side 
operators of e^'s form an inference chain ZI— XZ-ZX — 
IX - XY -YZ -YY - ZI. Note that this inference 
chain exhausts all the combination of operators on the 
right-hand side of eq.(Cl). Therefore, {e^ji is eventually 
periodic with period 7, or is eventually //. li 7 ] L, 
we must have = // for all i. We showed that the 
coefRcients of the operators touching / in O are all zero. 
We can repeat the argument to infer that O docs not 
involve stabilizer generators lying in between z — ±1. 

Consider the set of sites on the line given by a; = 
0,z = —1. A unit edge e[ in I' connecting (0, i,— 1) and 
(0,i-|-l,— l)isa linear combination of 



ZX -XZ 
XX - ZI 



XY - ZY 
ZZ - II. 



The inference chain is ZI-IY-XZ-YZ-YX- 
ZY — XX — ZI, which is again of eventual period 7. 
Therefore, O is a zero combination. The direction of the 
tunnel T can be along any coordinate axis, as is implied 
by three-fold symmetry. 

In a general case where 7 may divide L, we consider 
three tunnels Tx,Ty, intersecting at one 1x2x2 box 
B. It is easy to see that stabilizer generators {S) lying 
outside TxUTyUTz are independent; there are 2LxLyLz — 
8Lx — ALy — ALz + 16 independent generators. We can 
add to 5* more independent generators. Let the box B 
be given by < a; < 1, — 1 < y, ^ < 1. Consider a 
long contractible tube T' given by 1 < a; < L^r, < 
y < 1, — l<z<l. Then S' , the union of S and the 
set of stabilizer generators lying in the tube T', is a set 
of independent generators. Since T' contains A{Lx — 1) 
generators, \S'\ — 2LxLyLz—ALx — ALy — ALz + 12, which 
proves the upper bound k < 12L — 12. 



3. String logical operators 

There exist string logical operators for Code 0. By 
three-fold symmetry, these string logical operators form 
a triple. Using the notation in Table [TV] they are 



9^ - ZZ[1, 0, -l]p+iXI[l, 0, -l]pZZ[l, 0, -l]p_,, 
0^ ^ZX[-l, 1, 0]p+iXI[-l, 1, 0]pZX[-l, 1, 0]p_£, 
0^ ^ZY[0, -1, l]p+yXI[0, -1, l]pZY[0, -1, l]p^y, 

where z,x,y are directional unit vectors. We call them 
basic string logical operators for Code 0, or basic strings 
for short. Applying the techniques of deforming logi- 
cal string segments, one can show that any logical string 
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segment is a union of 'flat' logical string segment and 
some contiguous part of basic strings. Numerical result 
suggests that any long flat logical string segment (i^^ is 
trivial; it is trivial if the length of is greater than 3w, 
where w is the width of d°\ for w = 2, . . . , 600. There- 
fore, it is legitimate to conjecture that any long logical 
string segment is some product of basic strings. If this is 
true, we can consider a subsystem code by gauging out 
the logical qubits that are affected by basic strings. 

Interestingly, we can show that there are at least one 
logical qubit left in the resulting subsystem code in any 
finite periodic lattice Zj^. To see this, it is enough to 
calculate the commutation relation of independent ba- 
sic strings (up to the stabilizer group). After Gram- 
Schmidt procedure with respect to commutation sym- 
plectic form applied to independent set of basic strings, 
suppose we get 2fc;j-dimcnsional hyperbolic space, and fcj- 
dimensional isotropic space. Then the number of gauge 
qubits is kh + ki. 

is m 



It is not hard to see that 6*^ q 



(0,0,0) 



or 

(0,0,0) 



the stabilizer group 5o of Code in 1?^ ; it is equal to the 
product of QojQq inside the triangle formed by three 
basic strings. Also, 



z 

(0,0,0) 


f)Z 

(a,b,c) 


e5o. 


(0,0,0) 


f)X 




'(0,0,0) 


^{a.b,c 





for any a,b,c e such that a + b + c — 0. Therefore 
the maximally independent set of basic strings up to Sq 
is contained in 

r = {^^,j,o)'^{o,o,») ■ * ^ ^l}, 
i.e., 2kh + h< |r| = 2L. If L is odd, 

Joio] 



i 



)Z 

(0,i,0) 

)X — ^ 
(0,0, i) = ^ 



XI ^ ^0, 

[001] ^ <^ 
XI ^ '^0, 



where = means equality up to So- Hence 2kh+ki < 2L—2 

always commute, 

(0,i,0) 



for odd L. Note that 6*^ and 6*^ , 
and so do and ^fo.^'fi)■ 



Two basic strings 6*? 



and J Q-j anti-commute if and only if they meet at one 
site. We can write the commutation relation matrix of T 

^0 .'^ 




(over the binary field) as a;(r) 



where 



[i], [i±2] are the equivalence classes of integers modulo L, 
and 5 is the Kronecker delta. Note that rank(w) = 2fc/i. 
Each row of w' is a translation of another, and containes 
exactly two I's separated by 4 (L > 4). 

If L is odd, then L and 4 are relatively prime, and we 
see that rank(ti;') = L—1. Since 2kh + ki < 2L—2, ki = 



and the number of gauge qubits is L — 1. If L is even, 
we need to distinguish two cases. When 4 | L, only L — 4 
rows are independent. It is verified easily if one cyclically 
rotate the rows of oj' such that = 1. When 4 -f L, there 
is a row that contains two I's that are 2 columns apart, 
and hence L — 2 rows become independent. In short, 
kh = rankw' = 4[L/4] -4, and therefore k^ + ki < L + A. 
Since k > L + 6q2{L), we conclude that there are at least 
1 qubit left after gauging out the basic strings. 



Appendix D: Derivation of 2D toric code 

We apply our construction of iion-CSS cubic codes to 
2D square lattice. We argue that the 2D toric code is the 
unique stabilizer code under the construction. A square 
stabilizer generator s has 4 corner operators, a; is a 4 x 4 
matrix. Imposing the condition that s's define a stabilizer 
code, we get the general form: 

/O i j 0\ 

i j 



j i 
\0 J I 0/ 



where i,j £ F2. In order to ensure a single site operator is 
the identity, we require that rank(w) — 2m. If i = j = 0, 
rank(a;) =0. If i = j = 1, then rank(cj) = 2 and the 
realization of cj is 



X 



Z 



X 



This is the generator for 45°-rotated 2D toric code. If 
i = l,j = 0, then rank(tj) — 4 and a realization of uj is 



XI 




IX 



ZI 



IZ 



XI 



If i = 0, j = 1, we get the same realization of uj that is 
90°-rotatcd. Note that first qubits in the odd numbered 
rows interact only with second qubits in the even num- 
bered rows, and vice versa. Thus, the code is equal to 
the non-interacting two copies of 2D toric code. If the 
periodic lattice has odd linear size, it is a doubly-folded 
toric code. 



Appendix E: Numerical methods 

We describe our algorithm calculating k. As we have 
seen from proofs of the upper bounds for k ( Table [TTT| ), we 
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need to know the number of independent stabilizer gen- 
erators in a given finite periodic lattice. Since there are 
two qubits per site, a Pauli operator on Z/^ is expressed 
as a 4L'^- component binary vector. Given that the Pauli 
operator is Z- or X-type, the number of components is 
divided by two. Since there are 2L^ stabilizer generators, 
we need to calculate the rank of 21? x binary matrix 
[7. The Gauss elimination is an efficient algorithm, but 
the matrix is quite large. A naive sparse matrix method 
would not be of much help because U may get a large 
number of non-zero components as Gauss elimination al- 



gorithm proceeds. 

The proof of the upper bound gives a natural order 
of the sites and generators such that a large sparse sub- 
matrix of [/ is in a row echelon form; order the sites and 
generators such that the independence of generators used 
in the proof of the upper bound is evident from the form 
of U . In this way, we can represent V using memory 
size O(L^). It is also advantageous in view of time com- 
plexity. Note that since the field is binary, there is no 
multiplication. The naive Gauss elimination would re- 
quire 0(1?) additions, while the better representation of 
U needs only 0{L^) additions. 
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